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ABSTRACT 



We have obtained new images and high-resolution (R ~ 22400) near-infrared 
(1.2400-1.2575 /im) spectra of each component of the brown dwarf binary GJ 569 Bab 
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using the Adaptive Optics facility of the Keck II telescope and the NIRSPEC spectrom- 
eter. These data have allowed us to improve the determination of the astrometric orbit 
and to measure radial velocities of the components. We have used the astrometric and 
spectroscopic measurements to derive the dynamical mass of each brown dwarf and the 
systemic velocity of the pair by means of a \ 2 fitting technique. From various consid- 
erations, the mass of each component is likely in the range 0. 034-0. 07OM (GJ 569 Bb) 
and 0. 055-0. 087M© (GJ569Ba). This implies that the mass ratio, q, of the binary 
is greater than 0.4, being the most likely value q = 0.75-0.85. Adopting O.O72M0 as 
the most conservative location of the substellar limit for solar metallicity, our analysis 
confirms GJ569Bb as the first genuine brown dwarf known without any theoretical 
assumption. We have compared the dynamical masses of GJ 569 Ba and Bb, and their 
effective temperatures and luminosities, to the predictions of state-of-the-art theoretical 
evolutionary isochrones, finding that models exhibit good performance in the regime of 
high substellar masses if the binary is about a few hundred million years old. However, 
the surface gravities of GJ 569 Ba (M8.5V) and Bb (M9V) derived from our spectral 
analysis (the observed data have been compared to the latest synthetic spectra) appear 
to be smaller than the values provided by the evolutionary models. 

Subject headings: stars: low-mass, brown dwarfs — stars: individual (GJ 569 B) - 
stars: fundamental parameters — binaries: visual 

1. Introduction 

Since the discovery of the first brown dwarfs (Nakajima et al. (1995); Rebolo, Zapatero Osorio, 
& Martin (1995)), with masses below the hydrogen burning mass limit (Kumar (1963); Grossman, 
Hays, & Graboske (1974); Burrows et al. (1993); D'Antona & Mazzitelli (1994)), many efforts 
have been devoted to study their physical properties and formation processes. The frequency 
and characteristics of binary and multiple brown dwarf systems and brown dwarf-star pairs may 
provide key information on stellar /substellar /planetary formation. With the advent of new ground- 
based observational techniques, like the Adaptive Optics and coronography, and the high spatial 
resolution imaging capability of the Hubble Space Telescope, it has been possible to resolve very low 
mass companions around stars and binary brown dwarfs (e.g., Martin, Brandner, & Basri (1999); 
Close et al. (2002a), (2002b); Potter et al. (2002); Burgasser et al. (2003); Gizis et al. (2003); 
McCaughrean et al. (2004)). One spectroscopic binary brown dwarf has also been discovered in 
the young Pleiades open cluster (Basri & Martin (1999)). All these surveys indicate that substellar 
binaries are not rare, and that brown dwarfs are commonly found in orbit around low mass stars. 
These are rather nearby systems; the visual binary brown dwarfs known so far typically present 
projected separations between 0.9 and ~10AU (orbital periods of 3-30 yr or larger), permitting 
orbital mapping over the coming decades. 
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Eclipsing binaries have often served as valuable tools for the validation of structure and evo- 
lutionary models (e.g., Delfosse et al. (2000)). Accurate substellar mass, radius and luminosity 
measurements provide a crucial test of our understanding of substellar physics. To the best of our 
knowledge, no eclipsing brown dwarfs are identified so far. In the absence of eclipsing substellar 
systems, nearby visual brown dwarf pairs become benchmark objects to evaluate theoretical models 
in the substellar mass regime, because astrometric and spectroscopic studies will yield highly ac- 
curate dynamical masses of the components, although the majority of these binaries tends to have 
long orbital periods. The binary GJ 569 Bab is one of the very low mass pairs with the shortest 
periods known to date (2.4 yr, Lane et al. (2001a)). It can be used to calibrate state-of-the-art 
evolutionary isochrones. 

GJ 569 Bab is a binary dwarf moving around the chromospherically active M2.5V-type star 
GJ569 A (Hipparcos distance of 9.8 pc, Perryman et al. (1997)). The presence of a "single", very 
red, proper motion companion at a separation of 5" from the primary star was first reported by 
Forrest, Skrutskie, & Shure (1988). Follow-up low-resolution spectroscopy obtained by Henry & 
Kirkpatrick (1990) confirmed its cool nature (it was classified as an M8.5V dwarf). Martin et al. 
(2000) resolved the companion into two separate objects (GJ 569 Ba and GJ 569 Bb) using Adaptive 
Optics (AO) observations with the Keck II telescope. More recently, we obtained AO images and 
low-resolution spectra of each of the components (Lane et al. (2001a), hereafter Paper I), which 
have allowed us to determine the astrometric orbital parameters of the binary (see Table 1) and the 
spectral types of GJ569Ba (M8.5V) and GJ569Bb (M9V). Similar astrometric results obtained 
from speckle interferometry were also reported by Kenworthy et al. (2001). In this paper, we 
present high-resolution spectra of the resolved components obtained at different orbital epochs and 
the derivation of the dynamical masses of GJ 569 Ba and Bb. For the first time, one brown dwarf is 
unambiguosly confirmed without any theoretical model. This has allowed us to carry out a critical 
evaluation of substellar evolutionary model predictions. 



2. Observations and data reduction 

2.1. Adaptive Optics imaging and improved astrometric orbit 

In addition to the astrometry presented in Paper I, we have observed the binary GJ 569 Bab on 
two more occasions on 2001 Jun 28 and Sep 01. We used the 10-m Keck II AO facility (Wizinowich 
et al. (1988)), and the slit- viewing camera (SCAM) associated with the NIRSPEC instrument 
(McLean et al. (1998)). The pixel scale of the 256 x 256 Rockwell HgCdTe array of SCAM was 
/ /0168. Exposures were taken in the K'-band. The observational strategy, data reduction, and the 
analysis of the astrometry are all described in Paper I. We list in Table 1 the most recent relative 
positions of GJ569Bb with respect GJ569Ba as a function of the Modified Julian Date (MJD). 

These new astrometric data along with the previous astrometry allow us to determine the orbit 
of GJ 569 Bab to a certain degree of accuracy by fitting a Keplerian model. A description of the 



-4 - 



procedure by which the best Keplerian solution is obtained is fully detailed in Paper I. The updated 
best-fit parameters and their la uncertainties are given in Table 2 (columns 2 and 3), where we also 
provide the previous measurements for comparison (columns 4 and 5). The agreement is excellent, 
and the new and "old" values are very similar within the error bars. We note that the uncertainties 
of the updated parameters are reduced as compared to Paper I. This suggests that the relative 
orbit of the GJ 569 Bab pair is now determined to a high degree of confidence. The visual fit orbit 
is shown together with all the available astrometry in Fig. 1. 

It is important to note that the astrometric solution of the orbit in Paper I has an ambiguity of 
180° in the determination of the position angle of the ascending node, fi. This parameter sets how 
the plane of the orbit points in or out of the plane of the sky. From purely positional measurements 
we are able to identify the two nodes of the orbit, but it is not known which is the ascending node. 
Because the argument of periapsis, uj, is measured from the position of the ascending node, it also 
shares a similar ambiguity. This uncertainty does not affect the determination of the other orbital 
parameters or the total mass of the binary, and is only removed with spectroscopic data. The radial 
velocities that we have measured for the GJ 569 Bab system (Section 3.1) confirm that the angle of 
the node published in Paper I corresponds to the position of the descending node. In column 2 of 
Table 2 we provide the correct mean values for Q and to. 

From the updated astrometric orbit of the pair, the best-fit values of the period and semi-major 
axis correspond to a total mass of O.125M0 with a 3a uncertainty of ±O.O15M0. This is a purely Ke- 
plerian dynamical result, and no evolutionary models are involved in the calculations. The 3a upper 
mass limit of GJ 569 Bab is O.14OM . As discussed in Paper I, the similarity of the photometric and 
spectroscopic properties between GJ 569 Ba and Bb suggests that the mass ratio, q, of the pair is 
close to equal. This sets the 3a upper mass limit of the smallest object to be 0.070M Q (O.O65M at 
la), and the lower mass limit of the largest object to be 0.055M Q (0.060Mq at la). Magazzu, 
Martin & Rebolo (1993) did not detect lithium in the optical spectrum of GJ 569 Bab. Their data 
are dominated by the more massive component, which is brighter. Theoretical models (e.g., Bur- 
rows et al. (1997); Baraffe et al. (1998)) predict that objects more massive than 0.055M Q have 
efficiently depleted their atmospheric lithium at ages older than a few hundred million years. In 
this sense, models and our observations appear to be in agreement. We will discuss lithium in each 
component of the binary in a forthcoming paper. To investigate further the reliability of the mod- 
els, we must know the individual masses of GJ 569 Ba and Bb. The radial velocity measurements 
of Section 3.1 shall constrain better the mass of each component, providing a better estimate of q. 



2.2. Adaptive Optics high-resolution spectroscopy 

We have obtained J-band spectra with the Keck II telescope, the AO correcting system, and 
NIRSPEC, a cross-dispersed, cryogenic echelle spectrometer employing a 1024 x 1024 ALADDIN 
InSb array detector. These observations were carried out on four different occasions within the same 
orbital period during 2001; three of them coincide with astrometric observations. In the echelle 
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mode, we selected the NIRSPEC-3 (J) filter, and an entrance slit width of (X'43 (i.e, 3 pixels along 
the dispersion direction of the detector). The length of the slit was 12". For the present study 
we used only the spectral order that covers the following wavelength range: 1.2400-1.2575 /xm. 
This instrumental setup provided a nominal dispersion of 0.179 Apix" 1 (4.3 kmpix" 1 ), and a final 
resolution of 0.55 A, which corresponds to a resolving power i?<~ 22500 at 1.2485 fim. The journal 
of observations is given in Table 3 and shows the observing dates, exposure times and the air mass 
range of the observations. 

The slit was aligned with GJ 569 Ba and GJ 569 Bb so that the two objects were observed 
simultaneously. Spectra were collected at two different positions along the entrance slit. In order 
to remove atmospheric telluric absorptions, the near-infrared featureless stars HR5567 (A0V), 
HR5931 (A0V) and SAO 101508 (F8III) were observed very close in time and in air mass (within 
0.1 air masses). White-light spectra obtained with the same instrumental configuration were used 
for flat-fielding the target data. The dispersion solutions for the May, June and September data 
were derived from the arc lamp lines of Ar, Kr, and Xe, which were systematically taken after 
observing GJ 569. No calibrating arc lamps were observed in August. 

Raw spectra were reduced following conventional techniques in the near-infrared and using 
the REDSPEC software, which was particularly designed at the University of California for the 
analysis of NIRSPEC data. We also reduced and extracted spectra using the packages TWODSPEC 
and ONEDSPEC within the IRAF 1 environment, obtaining very similar results. Nodded images 
were subtracted to remove sky background and dark current. We selected different apertures to 
extract the individual spectra (between 4 and 6 pixels for GJ569Bb, and between 4 and 8 pixels 
for GJ569Ba). Figure 2 shows the cut along the spatial direction of one of our two-dimension 
spectra. As can be seen from the location and distribution of the apertures in the Figure, the 
cross-contamination between the extracted data of the two components is less than ~10% (except 
for the September data). We have checked that the relative radial velocity measurements obtained 
for the various apertures are consistent with each other within the error bars (see discussion of 
Sect. 3.1). Larger extraction apertures are useless for radial velocity measurements because the 
spectrum of each component of the binary is significantly contaminated by the presence of the 
companion. Extracted spectra were divided by normalized flat-fields, and calibrated in wavelength. 
The la dispersion of the fourth-order polynomial fit was 0.02 A (3.6% of the final resolution). 
We also removed some fringing that is internal to NIRSPEC and appears sometimes when using 
this instrument behind the AO system. The spectra of hot stars were used for division into the 
corresponding science spectra. These stars do not show any intrinsic atmospheric feature in the 
spectral coverage of the echelle order at 1.2485 fim. Finally, we multiplied the science spectra by 
the black body spectrum for the temperatures of 9480 K and 11200 K, which corresponds to the 
A0V and F8III classes, respectively (Allen (2000)). 



1 IRAF is distributed by National Optical Astronomy Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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The final spectra are shown in Fig. 3. No heliocentric or radial velocity corrections have been 
applied. The August data are not displayed because they lack a reliable wavelength calibration. 
The most noticeable spectral features are due to the absorption Ki resonance doublet at 1.2436 
and 1.2526 /im. The spectrum of the early M-type star GJ569A shows narrow atomic lines and 
no other relevant features are apparent. However, the Ki lines of the late-M dwarfs GJ569Ba and 
C appear quite broad, as expected for rather cool atmospheres. We note that these atomic lines 
are moderately stronger in G J 569 Bb than in Ba, which is indicative of the slightly cool nature of 
the smallest object. Many other absorption features are likely due to H2O, TiO, FeH and CrH, as 
discussed in Sect. 3.2. 

There are clear advantages in determining radial velocities from the spectra of Fig. 3. On the 
one hand, this wavelength interval is almost completely free of strong Earth's telluric absorption 
lines that might introduce uncertainties in the Doppler shift calculation. There are a few of these 
telluric lines in the red part of the spectra. Nevertheless, we are confident that the science spec- 
tra have good cancellation of these terrestrial atmospheric features. In addition, the wavelength 
solution of the spectral order at 1.2485 fj,m is reliable because it is based on seven arc lamp lines 
distributed all across the dispersion axis of the detector. Finally, the science spectra have relatively 
good signal-to- noise ratio at these wavelengths, and display numerous characteristic features that 
may yield accurate radial velocity measurements. 

3. Analysis 

3.1. Radial velocities 

The Doppler shifts of GJ 569 Bb were determined by a Fourier cross-correlation technique, 
using the spectrum of GJ 569 Ba as a template. The relative velocities of GJ 569 Bb cross-correlated 
against GJ 569 Ba obtained for the four observing nights are given in the second column of Table 4. 
According to the measurements, GJ 569 Bb is moving toward us, while GJ 569 Ba is moving further 
away at the time of the observations. This has allowed us to determine unambiguosly the position 
of the ascending node angle of the orbit of the pair (Table 2). We have measured the observed radial 
velocities of each component from the centroids of the K I lines, obtained by marking two continuum 
points around the lines and fitting the profiles within iraf. These velocities were corrected for the 
rotation of the Earth, the motion of the Earth about the Earth-Moon barycenter, and the orbit 
of the barycenter about the Sun to yield heliocentric radial velocities, which are listed in columns 
3 and 4 of Table 4. Heliocentric velocities cannot be derived from the August data because the 
spectra lack a reliable wavelength calibration. 

The precision of all the velocities may be judged by the standard deviation of the measurements 
for different regions of the spectra. The cross-correlation technique was able to achieve precisions of 
about 1/4 pixel (an accuracy of the order of ±1.1 kms -1 ) for good signal-to-noise data. The velocity 
accuracy became a factor of 1.5 to 2 worse when analizing the data of GJ569Bb, because this 
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object is fainter than GJ569Ba and its spectra show lower signal-to-noise ratios (Fig. 3). Similar 
uncertainties are affecting the radial velocities derived from the centroids of the atomic lines. We 
have measured the heliocentric velocity of the star GJ569A using its 2001 June spectrum, and 
have found it to be — 8.0 ±1.1 kms" 1 . Various authors have investigated the velocity variability of 
GJ 569 A during the last two decades, concluding that this is a rather stable star over a period of ten 
years with a mean heliocentric velocity of —7.17 ±0.28 km s -1 (Marcy, Lindsay & Wilson (1987); 
Marcy & Benitz (1989)), -8.6 ±1.0 km s" 1 (Barbier-Brossat & Figon (2000)), and -7.2 ±0.4 km s" 1 
(Nidever et al. (2002)). Our measurement and its error bar are in agreement with the values seen 
in the literature. We are confident that the velocities of Table 4 are accurate to the claimed 
uncertainties. 



3.2. Atmospheric gravity and rotation from spectral synthesis 

By comparing our observed spectra to synthetic energy distributions, we can obtain information 
on the atmospheric temperatures, gravities and rotation of GJ 569 Ba and Bb. Theoretical spectra 
were computed in a classical framework, i.e., assuming LTE, plane-parallel media, and no sinks 
and sources of energy in the atmosphere. We used the WITA6 code (Pavlenko (2000)), which 
takes into account a total of 100 species. The equations of ionisation-dissociation equilibrium were 
solved for a plasma made of neutral atoms, ions and molecules. The constants of the formulas for 
chemical balance were taken from Tsuji (1973) and Gurvitz et al. (1982). Line lists come from 
different sources, depending on the species. Atomic line lists were taken from VALD (Kupka et al. 
(1999)), lines of TiO from Plez (1998), CN from Kurucz (1993), CrH and FeH from Burrows et 
al. (2002) and Dulick et al. (2003), respectively, and H 2 from Partrige & Schwenke (1997). The 
partition function of H2O was computed as in Pavlenko (2002) to properly account for the splitting 
of the lines. We adopted the Voigt profile for the shape of each molecular and atomic line; the 
damping constants were computed using data from the original databases, and if no values were 
published, computed following an approximation of Unsold (1955). Due to the low temperatures of 
the atmospheres of late-M dwarfs, Stark broadening may be discarded, whereas pressure broadening 
prevails. 

A grid of 20 synthetic spectra was computed for solar metallicity, dusty model atmospheres 
of Allard et al. (2001), effective temperatures, T eff , of 2600, 2400, 2300, 2000 K, and gravities 
log g = 4.0, 4.5, 5.0, 5.5, 6.0 (cms~ 2 ). We have also computed grids for the COND model atmospheres 
of Allard et al. (2001), and the very recent ucm models of Tsuji (2002) for the critical temperature 
of 1800 K. The three sets of models differ in the treatment of dust (see the references above). Late- 
M and early-L dwarfs are supposed to have dusty upper atmospheres, with dust particles affecting 
the thermal structure. Hence, the COND models, which assume a complete dust-segregation, are 
not a good approximation for the case of GJ 569 Bab. Theoretical spectra were obtained in the 
wavelength range of our observations, with wavelength step of 0.5 A, and with microturbulent 
velocity of 2kms~ 1 . Element abundances are as in Anders & Grevesse (1989). The contributions 
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of various molecular and atomic species to the total opacity in the 1.2400-1.2575 /im spectral region 
are shown in Fig. 4, from which we see the absorption due to K I dominates. The remaining features 
of the observed spectra are likely due to H2O, CrH and FeH. However, the absorption at 1.2462 /im 
does not appear to have a counterpart in any of the molecular data of the figure. 

The effect of changing T e g- and log g in our computations is shown in Fig. 5. We note that 
CrH was excluded from our final synthetic spectra, because it appears too strong in the spectral 
region of the figure (see Fig. 4). Our spectroscopic data are better reproduced by computations 
that do not include the opacities of this molecule. As seen from Fig. 5, temperature has a relatively 
larger impact on the energy distributions than has gravity. It is also obvious that the Ki lines 
are remarkably sensitive to both temperature and gravity, i.e., the doublet becomes systematically 
weaker at high temperatures and low gravities (e.g., see McGovern et al. (2004)), while the main 
molecular species at these wavelengths do not appear to change notably with gravity. This justifies 
the use of the Ki atomic lines to estimate the surface gravity of GJ569Ba and Bb. 

Our observed spectra have been compared to the various grids of synthetic data, which were 
computed including the broadening of the alkali lines due to van der Waals collisions with H and 
H2 (important in cool, neutral atmospheres). The van der Waals treatment is an approximation 
applicable to the Ki doublet at 1.25 p, because these are weak lines (e.g., Burrows, Marley, Sz 
Sharp (2000)). The comparison has been performed as in Jones et al. (2002), Pavlenko & Jones 
(2002) and Pavlenko et al. (2004). The analysis procedure described by these authors is based on 
a minimization technique that provides the best-fit T e ff, logg and v TOt s'mi (line broadening due to 
atmospheric rotation), where i is the inclination angle of the rotation axis. Of the three sets of 
model atmospheres, the dusty models yield the smallest values of the minimization function. The 
best-fit parameters are as follows: for both GJ 569 Ba and Bb, T eS = 2400 K and log 5 = 4.5 (dusty 
models), and T e g- = 2500K and logg = 4.5 (UCM models). Tsuji's (2002) models provide higher T e g 
because they predict lower temperatures at the atmospheric layers where lines are produced. A 
plot depicting the (T c g, log g) contours of the minimization function is shown in Fig. 6, and Fig. 7 
displays the best-fit synthetical spectra together with the observed data. We note that while the 
alkali lines are nicely reproduced by the theoretical data, the molecular absorptions are not well 
matched by the computations. 

Unfortunately, we do not have a reference slow-rotator spectrum (of similar spectral type 
and observed with the same instrumentation) to measure the rotational velocities of G J 569 Ba 
and Bb. The fitting technique of Pavlenko et al. (2004) provides estimates for these velocities. 
Despite the unknown inclination angles of the two objects, our spectral analysis suggests that 
GJ569Ba and Bb are fast rotators, with Ba (v TOt suu ~ 37 ± 15 kms -1 ) possibly rotating slightly 
faster than Bb (v rot sinz ~ 30 ± 15 kms -1 ). We note that there is a large uncertainty associated 
to these measurements, mainly due to the poor knowledge of the molecular opacities and the 
uncertainty in the T e g- determination, which is about 100 K. A moderately cooler temperature than 
that of G J 569 Ba and a similarly high rotational velocity may explain the K 1 lines broadening 
and the dissimilarities of GJ569Bb. These rapid velocities are not surprising since many field 
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low-mass stars and brown dwarfs display high values of v ro t sin i as indicated by the spectroscopic 
measurements of Basri &: Marcy (1995), Tinney & Reid (1998), Basri (2001) and Reid et al. 
(2002). According to evolutionary models (e.g., Baraffe et al. (2003)), the radius expected for 
0.060-0.070M Q brown dwarfs is in the range 0.130-0.092 R Q for ages between 300 and 1000 Myr. 
This would imply short rotation periods between 3sini and 5sini hours. Some field and cluster 
brown dwarfs of similar spectral classes have been identified with such rapid rotations, like Kelu 1 
(1.8 h, Clarke, Tinney, & Covey (2002)), and SOri27, SOri28 and SOri45 of the a Orionis cluster 
(2.8, 3.3 and ~3h, respectively, Caballero et al. (2004); Zapatero Osorio et al. (2003)). 

Because the rotation velocities of GJ 569 Ba and Bb appear to be different, and assuming that 
both objects might have similar inclination rotation angles, it is unlikely that the two components 
are tidally pseudo-synchronized (i.e., the angular velocity of rotation is similar to the angular 
velocity of the companion crossing the sky during periastron passage). This contrasts with the 
case of the Pleiades member PP115 (120 Myr), which is, to the best of our knowledge, so far the 
only known spectroscopic binary brown dwarf (Basri k, Martin (1999)). The rotation of each of 
the components is similar and slow (lOkms -1 ), suggesting that the two objects are tidally pseudo- 
synchronized regardless their youth. The very short (the semi-major axis is 30 times smaller than 
that of GJ 569 Bab) and eccentric (e = 0.42) orbit of PP115 may also contribute to its pseudo- 
synchronism. 

We have also compared our spectroscopic data of GJ 569 Ba and Bb to computations carried 
out considering line broadening due to collisions with H (collisions with H2 were excluded). This 
comparison yields higher values of the minimization function, the same T c g, similar rotational 
velocities within 5 km s _1 , and higher log g values (log g = 5.0), for both the dusty and UCM models. 
Hence, we will adopt log g = 4.5-5.0 as the best surface gravity estimates for GJ569Ba and Bb. 
This result is in rough agreement with the recent calculations done by Gorlova et al. (2003). 
Concerning T e Q, the range 2400-2500 K is indeed very close to the estimates of Paper I, where 
we found 2440 ± 100 K for GJ569Ba and 2305 ± 100 K for GJ569Bb. We will use these values 
throughout the present paper. 

4. The mass of GJ 569 Ba and G J 569 Bb 

Radial velocities combined with the knowledge of the orbital parameters from astrometry 
robustly constrain the individual masses of visual binary members. However, in the case of the pair 
GJ 569 Bab there is an additional parameter that has to be considered: its orbital motion around 
the center of mass of the triple system GJ 569 ABab. 

The relative velocity of GJ 569 Bb with respect to Ba (or vice versa) is independent of any other 
velocity. It is only a function of the orbital parameters and the total mass of the pair. We have 
plotted in Fig. 8 our velocities (solid circles) against the epoch of the observations. Overplotted 
is the velocity curve computed for the updated orbit of Table 2. The agreement between the 
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spectroscopic and astrometric data is excellent (within la the observational error bars), confirming 
the good quality of the Keplerian solution. 

The mass of each member of the binary is related to the radial velocity by the following 
equations: 

K = 2irasmi 

P(l- e 2)i/2[M Bo (l + <?)] { ' 

VBa = v + K q Mbh [cos [y + lo) + ecosu] (2) 

VBb = Vo - K M Ba [cos (v + uj) + ecoso;] (3) 

where P, i, a, w, e are the orbital parameters as given in Table 2, v is the true anomaly, and the 
velocity and mass of GJ569Ba are indicated by vs a and Me a , respectively. Equations 1 to 3 are 
formulated using the mass ratio, q = Msb/Msa, of the binary. The quantity v Q accounts for the 
long-term value of the radial velocity of the centre of mass of the binary, which includes the radial 
velocity of the triple system, the orbital motion of the centre of mass of the pair around GJ 569 A, 
and the average velocity of the binary about its barycentre. The total mass of GJ 569 Bab and the 
shape of the velocity curves of each component are fixed by the astrometric orbit; hence, there are 
only two variables in equations (2) and (3): Ms a and v . 

We have produced a grid of pairs of computed velocity curves for Ms a and v values in the 
range 0.055-0. lOOM© and from —9 to — 15kms _1 , respectively, and have compared these velocities 
to our observed heliocentric data of Table 4. We note that O.O55M is the minimum possible mass 
of GJ569Ba (at 3a) according to purely astrometric and photometric considerations (Section 2.1). 
The upper limit to the mass has been arbitrarily imposed. The trial long-term velocity interval 
covers the range between the largest and smallest heliocentric velocity of the components. To find 
the optimal mass and long-term velocity we have minimized the chi-squared (x 2 ) function given by 
the following equation: 



v - v. 



12 



X J = S^-^ (4) 



a, 



where ai stands for the errors in the velocity measurements. This is, observed velocities (v{) 
were weighted by their corresponding uncertainties so that the best quality measurements are more 
important in finding the best fit. Unweighted heliocentric velocities were also compared to the com- 
putations by imposing o"j = lkms~ 1 . Figure 9 shows the resulting x 2 against mass, and Table 5 
summarizes the best fits. The best-fit radial velocity curves correspond to the most likely masses 
of 0.071M© (GJ569Ba), O.O54M (GJ569Bb) for the weighted data, and O.O68M (GJ569Ba), 
O.O57M (GJ569Bb) for the unweighted data. We note, however, that six radial velocity measure- 
ments (three per object) are probably too few to accurately constrain the systemic velocity of the 
binary. Consequently, the uncertainty of the parameters including the covariance of v and Mj^ a is 
rather large, being ±1.8 kms -1 and ±O.O52M , respectively, obtained by increasing the \ 2 by 2.3. 
Dispersions of ±0.45 kms -1 and ±0.011M Q (given in Table 5) are derived from the various x 2 solu- 
tions found when changing the astrometric orbital parameters by their corresponding uncertainties 
(given in Table 2). This mass uncertainty is also consistent with the discussion below. 
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There is another way to constrain the mass ratio of the binary. It makes use of the luminosity 
and T e g of each component. We will update the luminosities of Paper I by using the most recent 
J-band bolometric corrections given by Dahn et al. (2002). GJ569Ba and Bb were assigned 
updated luminosities of logL/L = —3.35 ± 0.07 and — 3.58±0.07, respectively. By applying the 
astrophysical definition of effective temperature, we have obtained the radii of the components and 
found them to be: R Ba = 0.119 ± 0.020 R Q and R B b = 0.102 ± 0.020 R Q . The uncertainty in the 
radii has been derived from the errors of ±0.07 dex and ±100 K in luminosity and T e g, respectively. 
In Sect. 3.2 we have proved that the surface gravities of GJ 569 Ba and Bb are quite alike, implying 
that the mass ratio of the pair, q, is directly related to the square ratio of the radii. Thus, it 
follows that q = 0.74 ± 0.25. This result is indeed very similar to that obtained from the radial 
velocity fitting (Table 5), suggesting that the long-term velocity of GJ 569 Bab and the masses of 
the components are reliable to some extent. 

As discussed in Sect. 2.1, the 3 a upper limit of the mass of GJ569Bb is O.O7OM . At the 
99% confidence level, the dynamical mass of GJ569Bb is in the interval 0. 034-0. 7OM , and the 
dynamical mass of GJ569Ba is in the range 0. 055-0. 087M . These intervals nicely overlap with 
the mass estimates listed in Table 5 of Paper I, which were based solely on the comparison of the 
photometry with evolutionary models. More radial velocity measurements are needed to decrease 
the mass uncertainty and detect the presence of much tiny objects orbiting any of the brown dwarfs. 
This sets the mass ratio of the binary, q, to be >0.4, being the most likely value q = 0.75-0.85. The 
most conservative location of the hydrogen burning mass limit, i.e., the borderline between stars 
and brown dwarfs, is O.O72M for solar metallicity (Chabrier & Baraffe (2000)). The dynamical 
mass of GJ 569 Bb is definetively below the substellar frontier, confirming it as the first genuine 
brown dwarf known independently of models. GJ 569 Ba is very likely a brown dwarf too, although 
it may lie at the star-brown dwarf borderline. 

5. Discussion 

For the first time, it is possible to evaluate the predictions of theoretical evolutionary models 
in the substellar regime, by comparing isochrones to the dynamical mass, T e g, gravity and lumi- 
nosity of GJ569Ba and Bb. For this independent test, we need to constrain the age of the triple 
system. On the other hand, GJ 569 ABab may provide a calibration of the brown dwarf-mass-age- 
luminosity relation. 

5.1. The age of GJ 569 ABab 

The membership of GJ 569 ABab to a stellar kinematic group would help us set limits upon the 
age of the system. The space velocity components, UVW, of GJ569A put it firmly in the young 
disk region of the space motion groups as defined by Eggen (1969). Furthermore, these components 
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can be compared to the space motion of the Ursa Major (UMa) moving group (core and extended 
halo, sometimes called the Sirius supercluster) and the new supercluster found by Chereul, Creeze, 
& Bienayme (1999), which has a mean velocity between the Sun and the Sirius supercluster. 

In a recent paper, King et al. (2003) used the radial velocity of GJ 569 A (—7.17 ± 0.28 kms -1 ) 
measured by Marcy et al. (1987) and the Hipparcos parallax and proper motion of GJ569A to 
study its Galactic velocity components. These authors concluded that the membership of GJ569 A 
in the young UMa moving group is uncertain. Using the same data, Kenworthy et al. (2001) argued 
that the kinetics of GJ569 A is within 7.2 km s -1 of the space motion of UMa. Nevertheless, none 
of these authors took into account the orbital motion of G J 569 A about the barycenter of the triple 
system. The mass of this M2.5V-class star is estimated at 0.3-0.4M Q (Gorlova et al. (2003)). By 
applying simple Keplerian equations, we inferred that the observed radial velocity of GJ 569 A may 
be affected by an uncertainty of up to ±1.3kms _1 . Hence, the velocity vector of GJ569A and its 
1 a error bar turn out to be U = 7.8 ± 0.7, V = 3.3 ± 0.3, W = -13.3 ± 1.1 kms -1 . We have used the 
mathematical formulas of Johnson & Soderblom (1987) to calculate the space- velocity components 
and their uncertainties from the errors in the observational quantities (parallax, proper motion and 
radial velocity). Our result is in very much agreement with that of Kenworthy et al. (2001) and 
King et al. (2003), except for our larger error bars that mainly come from the uncertainty in the 
orbital velocity of GJ 569 A. 

From the literature, we have obtained that the mean Galactic space-velocity components of 
the UMa moving group are the following: [7 = 12.8, V = 2.2, W= -8.6kms _1 (Eggen (1992); 
Soderblom k Mayor (1993); Orlov et al. (1995); Chen et al. (1997); Asiain et al. (1999); Chereul 
et al. (1999); Montes et al. (2001); King et al. (2003)). The velocity dispersion of UVW among 
individual stars of the UMa group is typically around 3kms _1 (Madsen, Dravins, Sz Lindegren 
(2004)), although dispersions as high as 7kms~ 1 are seen in various works. This large scatter is 
not surprising since the group is unbound and spatially extended in the Galactic disk. The space 
motion vector of G J 569 A overlaps with that of the UMa group if the uncertainties are multiplied 
by a factor of 2, which indicates that its membership in the UMa cluster is possible, but not very 
probable. In the kinematic planes shown in Fig. 17 of Chereul et al. (1999), GJ569 A is positioned 
to the left of the Sirius supercluster, very close to the location of the new supercluster reported by 
these authors. The UVW coordinates of GJ 569 A differ at only the 1 a level with respect the space 
motion of this new supercluster, in particular with the stream moving at U = 5.3 ± 2.7, V = 6.5 ± 3.6, 
W = — 11.4 ± 4.0 kms -1 (stream 2-43 according to the nomenclature of the discoverers). On the 
basis of its space velocity, GJ 569 may kinematically belong to any of these two moving groups. 

The age of the UMa cluster has been widely discussed in the literature. The canonical quoted 
age of the nucleus is around 300 Myr (Soderblom et al. (1993), and references therein), although 
some older values (up to 600 Myr) are claimed more recently (King et al. (2003), and references 
therein). The massive stars of the new supercluster of Chereul et al. (1999) have ages very similar to 
the Sirius supercluster (100-1000 Myr). Nevertheless, based on Stromgren photometry, the stream 
2-43 contains rather young stars, with ages in the interval 300-800 Myr. Hence, we conclude that 
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if GJ 569 is a member of any of these kinematic groups, its age must be constrained in the range 
300-800 Myr. 

The spectroscopic and photometric properties of the star GJ 569 A also support the youth of 
the system. The age of GJ569A has been deeply studied in the literature (e.g., Forrest et al. 
(1988); Henry & Kirkpatrick (1990); Martin et al. (2000); Kenworthy et al. (2001); Gizis, Reid, 
& Hawley (2002)). The various groups conclude that GJ569A exhibits many of the attributes 
of young stars (chromospheric activity, flaring), and claim that the likely age of GJ569A is in 
the interval 120-1000 Myr. In addition, the color-magnitude diagrams shown in Paper I, which 
included the decomposition of the near-infrared photometry of the binary, also suggest that GJ 569 
is a young system. 

5.2. Mass luminosity and mass -T c r relations 

We note that our T e g- values from Paper I are within 50 K of those recently found by Gorlova 
et al. (2003), and within 75 K of the new T e ff-spectral type calibration by Dahn et al. (2002). The 
updated luminosities of Sect. 4 are also within the error bars of previous estimates. 

Figures 11 and 12 show GJ569Ba and Bb plotted on the mass-luminosity and mass-T e g 
diagrams, respectively. Both the weighted and unweighted dynamical masses are displayed. For 
the figures, we have considered four different sets of pre-main sequence evolutionary models: the 
NextGen models of Baraffe et al. (1998), the corrected isochrones of D'Antona & Mazzitelli (1997), 
hereafter DM97, the "Arizona" tracks of Burrows et al. (1997), and the most recent release of 
the COND models by Baraffe et al. (2003). These include essentially all of the models widely 
used by many authors for substellar studies. The two evolutionary tracks of Baraffe et al. (1998, 
2003) are usually referred to as the "Lyon" models. For a detailed description of the evolutionary 
models, see the corresponding authors. Nevertheless, we note that all of them are state-of-the- 
art calculations that employ sophisticated physiscs and boundary conditions. Since the two brown 
dwarfs of GJ 569 Bab are supposed to be coeval, a single isochrone is expected to fit both components 
simultaneously. Models must also reproduce the absolute and relative location of the brown dwarfs 
in the diagrams. 

From Figs. 11 and 12, the four sets of isochrones indicate that GJ 569 Bab is relatively young, 
~200-1000 Myr, with a likely age at around 300 Myr. We note that this result agrees with the 
discussion on the system age of the previous section. It is important to remark that all models 
yield reasonable fits to the observed parameters of the brown dwarf binary only if the system turns 
out to be about a few hundred million years old. In the following, we will briefly discuss on the 
fine details of the comparisons. 

We recall that the masses of GJ 569 Ba and Bb have been derived independently of models, 
whereas temperatures and luminosities hinge to some extent on external calibrations. For the Lyon 
and Arizona models, both the primary and secondary components of the binary lie within 1 a of 
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the 300-Myr isochrones depicted in Fig. 11. We are now considering the unweighted dynamical 
masses. However, the luminosity prediction of the 300-Myr DM97 track for the dynamical mass of 
GJ 569 Bb deviates by 2 a. In general, the Lyon and Arizona theoretical mass-luminosity relations 
for ages around a few hundred million years agree acceptably well with our observations of the 
binary GJ 569 Bab at the level of 1 a the observational uncertainties. 

Larger differences are found when comparing G J 569 Ba and Bb to the mass-T e g relations 
shown in Fig. 12. While the Lyon models at the age of 300 Myr appear to reproduce within 1 a 
the slope formed by the members of the pair, GJ 569 Ba and Bb lie at different isochrones of the 
Arizona models. The secondary and primary appear between 1 and 2 a up and down, respectively, 
of the predictions of the 300 Myr DM97 isochrone. 

5.3. Model atmospheres and evolutionary isochrones 

For a test of consistency between evolutionary models, model atmospheres and T e s calibrations, 
we have obtained the surface gravity of GJ 569 Ba and Bb by means of synthetic spectral fitting 
(see Sect. 3.2). Radii can be derived from logg and the Newton's law for Gravitation. Adopting 
g Q = 27380.215 cm s~ 2 and the unweighted dynamical masses, the radius of GJ569Ba and Bb is 
R/Rq = 0.18 1q 05 an d 0-17 1qo4' respectively. These values are plotted against mass in Fig. 13, 
together with the Lyon and Arizona isochrones. There is a large discrepancy (about a factor of 
1.6) between the predictions of the evolutionary models and the radii obtained from the spectral 
fitting. Such a discrepancy translates into a difference in log g of the order of 0.5 dex. Very recently, 
Mohanty et al. (2004) have reported finding very low gravities in young late-type stars of the Upper 
Scorpius association as compared to theoretical evolutionary model predictions. According to these 
authors, the differences are the largest for the coolest types (>M7). 

We note, however, that the radii derived in Sect. 4 are in much better agreement with the 
evolutionary models if the age of the pair is about a few hundred million years. 

The discrepancy between the radii predictions for GJ 569 Ba and Bb and those inferred from 
the spectroscopically derived surface gravities could be due to a combination of the uncertainties 
in the T e g-spectral type calibrations, the spectral synthesis, the model atmospheres, the log g step 
employeed in the synthetic calculations, and the evolution isochrones (see Mohanty et al. (2004) 
for a rather extensive discussion). On the other hand, the formation processes and evolution of 
binaries are not well understood, and tracks describing the evolution of "single" objects might not 
be valid. Nevertheless, warmer values of the T c g assignement of GJ 569 Ba and Bb would bring to a 
better harmony the radii predictions of evolutionary models and spectral fitting techniques. From 
our spectral fitting analysis, we observed that larger best-fit values of logg are found for hotter 
temperatures (see the contours of Fig. 6), e.g., 2600 K and log # = 5.0 (collisions with H and H2) 
and 5.5 (collisions with H), although they do not minimize the function defined by Pavlenko et 
al. (2004). A temperature of 2600 K is marginally consistent with the most recent calibrations for 
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M8.5-M9 dwarf objects. On the other hand, an increase of 150 K in the T e g of GJ569Ba and Bb 
does not change significantly the discussion of Sect. 5.2. Actually, the Arizona models would fit 
better the observed parameters of the binary at 300 Myr, and the Lyon models would predict a 
slightly younger age for the system. 

In any case, the spectral fitting of the Ki doublet at 1.25 /im and the radii derivation of Sect. 4 
suggest that the radii of G J 569 Ba and Bb are larger than the evolutionary model predictions for the 
oldest ages (~1 Gyr). This might be due to the young age of the objects (the brown dwarfs would 
be contracting); on the other hand, part of this discrepancy might be real. Similar discrepancies 
have been observed in low-mass eclipsing binaries and in the interferometric radii measurements 
of M-dwarfs as reported by various groups. For example, models fail to predict the radius of the 
transiting planet HD 209458 b (predictions are too small, Baraffe et al. (2003)), and underestimate 
the size of cool stars up to 10% (e.g., Lane, Boden, & Kulkarni (2001b); Ribas (2003); Segransan 
et al. (2003)). For a definitive and critical constraint of substellar model predictions, the size of 
each brown dwarf of the binary GJ 569 Bab should be determined by means that do not include any 
theoretical calculations or external calibrations, e.g., near infrared interferometry, or photometric 
monitoring to determine rotation periods and infer radii by assuming that the rotation angles are 
perpendicular to the orbital plane of the binary. 

5.4. Final remarks 

The mass ratio of GJ 569 Bab is q > 0.4, very likely q = 0.75-0.85 (Table 5). Many other binary 
brown dwarfs have similar mass ratio estimates (based on the comparison to evolutionary models), 
like, e.g., PP115 (g = 0.87, Basri & Martin (1999)), 2MASS J0231016-040618 (g = 0.7, Close et 
al. (2002a)), 2MASS J1426316+155701 (q = 0.9, Close et al. (2002b)), and HD 130948BC (q = 0.9, 
Potter et al. (2002)). Close et al. (2003) presents the compilation of nearly all very low mass, 
field resolved binary systems known to date, concluding that they tend to have nearly equal mass 
components (^ ~ 0.9) . Furthermore, Martin et al. (2003) have investigated the binary fraction 
among young brown dwarfs of the Pleiades cluster (120 Myr), finding that the frequency of these 
binaries is about 15% (for separations between 7 and 12 AU), with a clear trend to high mass ratios 
(g>0.7). Nevertheless, the discovery of substellar pairs like elndiB (q = 0.6, McCaughrean et al. 
(2004)) suggests that binaries with moderate mass ratios are also possible, albeit they appear to be 
less frequent. Any model describing the formation mechanisms of very low mass stars and brown 
dwarfs should be capable of explaining these observations. 

6. Summary and conclusions 

Using the Keck II Adaptive Optics facility and the NIRSPEC instrument, we have obtained 
new near-infrared images that have allowed us to improve the astrometric solution of the substellar 
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binary GJ 569 Bab, which is a proper motion companion to the M2.5V-type star GJ569A. The 
improved orbital parameters lie within 1-a the previously published values of Lane et al. (2001a), 
suggesting that the orbit is now determined to a high degree of accuracy. The total mass of the 
pair turns out to be 0.125 ± 0.015Mq at the 3-a confidence level. 

We have also obtained high-resolution spectra of each component, GJ 569 Ba and Bb, using 
the echelle mode of NIRSPEC on four different occasions during one orbital period. The spectra, 
with a resolution of 0.55 A (R ~ 22500), cover the wavelength range between 1.2400 and 1.2575 fim, 
which includes the Ki doublet at 1.2436 and 1.2526 /im. We have used these spectroscopic data to 
measure radial velocities of each member of the binary system with an accuracy between 1.1 and 
2kms~ 1 . The combination of the astrometric orbit and the spectroscopic velocities has allowed us 
to derive the dynamical masses of G J 569 Ba and Bb and the systemic velocity of the pair by means 
of a x 2 fitting technique. We note that our results show a marked correlation, which yields relatively 
large uncertainties in the parameters. More radial velocities measurements are needed to break the 
covariance problem. Holding the long-term velocity of GJ 569 Bab fixed at its best-fit values, the 
masses obtained are Mqj 569 Ba = 0.068 and Mqj 569 Bb = O.O57M for the unweighted velocities, and 
Mgj 569 Ba = 0.071 and Mqj 569 Bb = 0.054Mq for the weighted data. This implies that the mass 
ratio of the binary is greater than 0.4, being the most likely value q = 0.75-0.85. The dynamical 
masses are affected by an uncertainty of 16-20% (the error increases if the covariance problem is 
taken into account). Because of this uncertainty, the more massive component, GJ569Ba may lie 
at the star-brown dwarf borderline. However, the mass of the secondary, GJ 569 Bb, is very likely 
below the substellar limit. For the first time, the substellar mass of one brown dwarf is dynamically 
confirmed, i.e., without any theoretical assumption. 

The photometric and spectroscopic properties of the primary star GJ569A, as well as its 
UVW space velocities, suggest that the age of the triple system is in the interval 300-800 Myr. The 
dynamical masses of GJ569Ba and Bb, along with their effective temperatures and luminosities, 
have been compared to the mass-luminosity and mass-T e ff relations given by four different sets of 
state-of-the-art evolutionary models. We have found that these models reasonably reproduce the 
absolute and relative locations of GJ569Ba and Bb in the luminosity vs. mass and T e g vs. mass 
diagrams if the age of the binary is about a few hundred million years. Hence, as for luminosities 
and Teg's, evolutionary models appear to exhibit good performance in the regime of high substellar 
masses. 

The observed high-resolution spectra of GJ 569 Ba (M8.5V) and Bb (M9V) have been compared 
to various spectral synthesis using three different sets of model atmospheres. Our spectral fitting 
procedure yields T e fj's that are very similar (within 100 K) to the calibrations for M8-M9-class 
objects available in the literature, and surface gravities that are quite small as compared to the 
predictions of evolutionary models. A similar discrepancy has been recently reported by Mohanty 
et al. (2004). We argue that slightly warmer values of the T e ff's of the binary components would 
bring to a better agreement the gravity predictions of evolutionary models and spectral synthesis. 



-17- 



We thank Barbara Schaefer, Randy Cambell and David Le Mignant for their help with the 
2001 Sep 01 Keck service observations. We also thank F. Allard and T. Tsuji for providing machine- 
readable versions of their model atmospheres. We thank the anonymous referee for his/her com- 
ments. Data presented herein were obtained at the W. M. Keck Observatory, which is operated as 
a scientific partnership between the California Institute of Technology, the University of California, 
and NASA. The Observatory was made possible by the generous financial support of the W. M. 
Keck Foundation. The authors wish to extend special thanks to those of Hawaiian ancestry on 
whose sacred mountain we are privileged to be guests. This research has made use of the Simbad 
database, operated at CDS, Strasbourg, France. YP acknowledges partial financial support from 
PPARC and the Royal Society and Small Research Grant from American Astronomical Society. 
BFL acknowledges support from a Pappalardo Fellowship in Physics. This work was partly carried 
out with support from the MCyT projects AYA2003-05355 and Ramon y Cajal. 

Facilities: Keck(NIRSPEC). 

REFERENCES 

Allard, F., Hauschildt, P.H., A., D.R., Tamanai, A., & Schweitzer, A. 2001, ApJ, 556, 357 

Allen, W. Bb. 2000, Allen's Astrophysical Quantities. Fourth edition, ed. Arthur N. Cox, New York: 
Springer Verlag, p. 151 

Anders, E., & Grevesse, N. 1989, GeGoAA, 53, 197 

Asiain, R., Figueras, F., Torra, J., & Chen, B. 1999, A&A, 341, 427 

Baraffe, I., Chabrier, G., Allard, F., k Hauschildt, P. H. 1998, A&A, 337, 403 

Baraffe, I., Chabrier, G., Barman, T.S., Allard, F., k Hauschildt, P.H. 2003, A&A, 402, 701 

Barbier-Brossat, M., & Figon, P. 2000, A&AS, 142, 217 

Basri, G., & Marcy, G. W. 1995, AJ, 109, 762 

Basri, G., & Martin, E.L. 1999, AJ, 118, 2460 

Basri, G. 2001, ARA&A, 38, 485 

Burgasser, A. J., Kirkpatrick, J.D., Reid, I.N., Brown, M.E., Miskey, C.L., & Gizis, J.E. 2003, 
ApJ, 586, 512 

Burrows, A., Hubbard, W.B., Saumon, D., & Lunine, J.I. 1993, ApJ, 406, 158 
Burrows, A., et al. 1997, ApJ, 491, 856 

Burrows, A., Marley, M., & Sharp, CM. 2000, ApJ, 531, 438 



-18- 



Burrows, A., Ram, R.S., Bernath, P., Sharp, CM., & Milsom, J. A. 2002, ApJ, 577, 986 

Caballero, J. A., Bejar, V. J. S., Rebolo, R., & Zapatero Osorio, M. R. 2004, A&A, submitted 

Chabrier, G., & Baraffe, I. 2000, ARA&A, 38, 337 

Chen, B., Asiain, R., Figueras, F., & Torra, J. 1997, A&A, 318, 29 

Chereul, E., Creze, M., & Bienayme, O. 1999, A&AS, 135, 5 

Clarke, F.J., Tinney, C.G., & Covey, K.R. 2002, MNRAS, 332, 361 

Close, L.M., Potter, D., Brandner, W., Lloyd-Hart, M., Liebert, J., Burrows, A., & Siegler, N. 
2002b, ApJ, 566, 1095 

Close, L.M., Siegler, N., Potter, D., Brandner, W., & Liebert, J. 2002a, ApJ, 567, L53 

Close, L.M., Siegler, N., Freed, M., & Biller, B. 2003, ApJ, 587, 407 

Dahn, C. C, et al. 2002, AJ, 124, 1170 

D'Antona, F., & Mazzitelli, I. 1994, ApJS, 90, 467 

D'Antona, F., & Mazzitelli, I. 1997, Mem. S. A. It., 68, 807 

Delfosse, X., Forveille, T., Segransan, D., Beuzit, J.-L., Udry, S., Perrier, C, & Mayor, M. 2000, 
A&A, 364, 217 

Dulick, M., Bauschlicher, C.W. Jr., Burrows, A., Sharp, CM., Ram, R. S., & Bernath, P. 2003, 
ApJ, 594, 651 

Eggen, O.J. 1969, ApJ, 155, 701 

Eggen, O.J. 1992, AJ, 104, 1493 

Forrest, W. J., Skrutskie, M.F., & Shure, M. 1988, ApJ, 330, L119 
Gizis, J.E., Reid, I.N., & Hawley, S.L. 2002, AJ, 123, 3356 
Gizis, J.E., et al. 2003, AJ, 125, 3302 

Gorlova, N.I., Meyer, M.R., Rieke, G.H., & Liebert, J. 2003, ApJ, 593, 1074 
Grossman, A. S., Hays, D., & Graboske, H. C. 1974, A&A, 30, 95 

Gurvitz, L.V., Weitz, I.V., & Medvedev, V. A. 1982, Thermodynamic properties of individual 
substances. Moscow. Nauka. 

Henry, T. J., & Kirkpatrick, J. D. 1990, ApJ, 354, L29 



-19- 



Johnson, D.R.H., k Soderblom, D.R. 1987, AJ, 93, 864 

Jones, H.R.A., Pavlenko, Ya.V., Viti, S., k Tennyson, J. 2002, MNRAS, 330, 675 
Kenworthy, M., Hofmann, K.-H., Close, L., et al. 2001, ApJ, 554, L67 

King, J.R., Villarreal, A.R., Soderblom, D.R., Gulliver, A.F., k Adelman, S.J. 2003, AJ, 125, 
1980 

Kumar, H. 1963, ApJ, 137, 1121 

Kupka, F., Piskunov, N., Ryabchikova, T.A., Stempels, H.C., k Weiss, W. W. 1999, A&AS, 138, 
119 

Kurucz, R. L. 1993, CDROMs 1-22, Harvard-Smisthonian Observatory 
Lane, B.F., Boden, A.F., k Kulkarni, S.R. 2001b, ApJ, 551, L81 

Lane, B.F., Zapatero Osorio, M.R., Britton, M., Martin, E.L., k Kulkarni, S. 2001a, ApJ, 560, 
390, Paper I 

Madsen, S., Dravins, D., k Lindegren, L. 2004, A&A, in press 
Magazzu, A., Martin, E.L., k Rebolo, R. 1993, ApJ, 404, L17 
Marcy, G.W., k Benitz, K.J. 1989, ApJ, 344, 441 
Marcy, G.W., Lindsay, V., k Wilson, K. 1987, PASP, 99, 490 
Martin, E. L., Brandner, W., k Basri, G. 1999, Science, 283, 1718 

Martin, E.L., Koresko, CD., Kulkarni, S.R., Lane, B.F., k Wizinowich, P.L. 2000, ApJ, 529, L37 
Martin, E.L., Barrado y Navascues, D., Baraffe, L, Bouy, H., k Damn, S. 2003, ApJ, 594, 525 
McCaughrean, M. J., et al. 2004, A&A, 413, 1029 

McGovern, M. R., Kirkpatrick, J.D., McLean, I. S., Burgasser, A. J., Prato, L., k Lowrance, P.J. 
2004, ApJ, 600, 1020 

McLean, I.S., et al. 1998, SPIE, 3354, 566 

Mohanty, S., Basri, G., Jayawardhana, R., Allard, F., Hauschildt, P., k Ardila, D. 2004, ApJ, in 
press 

Montes, D., Lopez-Santiago, J., Galvez, M. C, Fernandez-Figueroa, M. J., de Castro, E., k Cornide, 
M. 2001, MNRAS, 328, 45 

Nakajima, T., Oppenheimer, B.R., Kulkarni, S.R., Golimowski, D.A., Matthews, K., k Durrance, 
S.T. 1995, Nature, 378, 463 



-20- 



Nidever, D.L., Marcy, G.W., Butler, R.P., Fischer, D.A., & Vogt, S.S. 2002, ApJS, 141, 503 

Orlov, V.V., Panchenko, I.E., Rastorguev, A.S., & Yatsevich, A.V. 1995, AZh, 72, 495 (English 
transl. Astron. Rep., 39, 437 

Partrige, H., & Schwenke, D.W. 1997, J. Chem. Phys., 106, 4618 

Pavlenko, Ya. V. 1997, Astron. Rep, 44, 219 

Pavlenko, Ya. V. 2002, Astron. Rep, 46, 567 

Pavlenko, Ya.V., & Jones, H. R. A. 2002, A&A, 396, 967 

Pavlenko, Ya.V., Geballe, T.R., Evans, A., Smalley, B., Eyres, S.P.S., Tyne, V.H., & Yakovina, 
L. A. 2004, A&A, in press 

Perryman, M.A. C, et al. 1997, A&A, 323, L49 

Plez, B. 1998, A&A, 337, 495 

Potter, D., Martin, E. L., Cushing, M. C, Baudoz, P., Brandner, W., Guyon, O., & Neuhauser, R. 
2002, ApJ, 567, L133 

Rebolo, R., Zapatero Osorio, M.R., & Martin, E.L. 1995, Nature, 377, 129 

Reid, I.N., Kirkpatrick, J.D., Liebert, J., Gizis, J.E., Dahn, C.Bb., & Monet, D. G. 2002, AJ, 124, 
519 

Ribas, I. 2003, A&A, 398, 239 

Segransan, D., Kervella, P., Forveille, T., & Queloz, D. 2003, A&A, 397, L5 

Soderblom, D.R., Stauffer, J.R., Hudon, J.D., & Jones, B.F. 1993, ApJS, 85, 315 

Soderblom, D.R., & Mayor, M. 1993, AJ, 105, 226 

Tinney, C.G., & Reid, I.N. 1998, MNRAS, 301, 1031 

Tomkin, J., & Popper, D.M. 1986, AJ, 91, 1428 

Tsuji, T. 1973, A&A, 23, 411 

Tsuji, T. 2002, ApJ, 575, 264 

Unsold, A. 1955, Physik der Sternatmospheren, 2nd ed., Springer. Berlin 
Wizinowich, P.L., et al. 1988, Proc. SPIE, 3353, 568 

Zapatero Osorio, M.R., Caballero, J. A., Bejar, V. J.S., & Rebolo, R. 2003, A&A, 408, 663 



This preprint was prepared with the AAS IATgjX macros v5.2. 



- 21 - 



0.12 



0.08 



Jun '01 



Aug '99/ Sep '01 May '01 



O 

o 
Q 



0.04 







Feb '00 



Jan '01 



-0.04 



-0.08 



Jun '00 



" Aug '00 



-0.12 1 1 1 1 1 1 1 

0.12 0.08 0.04 -0.04 -0.08 -0.12 

RA (asec) 

Fig. 1. — The new astrometric data (2001 Jun and Sep) are plotted together with the best-fit 
orbit (dashed ellipse) and the previous astrometry (Lane et al. (2001a)). Error-bar crosses denote 
measurements and circles indicate the predicted location on the orbit at the time of the observations. 
North is up, and east is to the left. 
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Fig. 2. — Cut along the spatial axis of the 2001 May two-dimension spectra of GJ 569 Ba (brightest 
peak) and Bb. A total of 6 columns are averaged. The smallest and largest apertures used to 
extract the individual spectra are shown to the top. 
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Fig. 3. — NIRSPEC spectra of the GJ 569 triple system. The spectrum of the primary star GJ 569 A 
was obtained in 2001 Jun 28. The strong absorption features are due to the atomic Ki doublet at 
1.2436 and 1.2526 (im. Data are displaced upwards by one for clarity. 
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Fig. 4. — Contribution of different molecular and atomic species to the opacity in the spectral 
region of our observations. 
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Fig. 5. — Dependence of theoretical spectra on T e g (left panel) and log g (right panel). CrH is not 
included in these computations, and van der Waals broadening due to collisions with H and H2 is 
considered. Spectra, computed using the dusty model atmospheres of Allard et al. (2001), have 
been rebinned to the resolution of our data and rotationally broadened by 30kms _1 . An offset of 
0.5 units is added to each spectrum for clarity. 
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Fig. 6. — The contours of the "S" function, defined by Pavlenko et al. (2004), against T e g and logg 
(cms" 2 ) for the dusty model atmospheres. The best-fit parameters (4.5, 2400 K) are given by the 
asterisk. 
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Fig. 7. — Spectral data of GJ 569 Ba and Bb (thick lines) and the best-fit dusty theoretical 
spectrum (2400 K, log g = 4.5, thin lines). Note that while the alkali lines are nicely reproduced by 
the theoretical data, the molecular absorptions are not well matched by the computations. 
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Fig. 8. — Relative velocity of GJ 569 Bb with respect GJ 569 Ba plotted as a function of the epoch 
of the observations. Overplotted onto the data is the predicted velocity curve for the known orbital 
parameters and total mass of the pair. 
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Fig. 9. — The \ 2 as a function of possible masses of GJ 569 Ba. For each mass, the systemic 
velocity of the pair has also been minimized. The unweighted (left panel) and weighted (right 
panel) heliocentric velocities of G J 569 Ba and Bb are best reproduced by the mass corresponding 
to the minimum of the \ 2 (vertical solid lines). The mass step in the calculations is 5xlO" 4 M . 
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Fig. 10. — Heliocentric velocities of GJ569Ba (circles) and Bb (triangles), together with the best- 
fit velocity curves (solid and dashed lines). The dotted, horizontal line indicates the centre of mass 
velocity of the pair. 
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Fig. 11. — Mass-luminosity relationships from various groups (NextGen models of Baraffe et al. 
(1998) — top left; COND models of Baraffe et al. (2003) — top right; D'Antona & Mazzitelli (1997) 
models — bottom left; Arizona models of Burrows et al. (1997) — bottom right). The weighted 
(open circles) and unweighted (solid circles) masses of GJ 569 Ba and Bb are also plotted together 
with 1 a error bars. The vertical dotted line represents the substellar limit. Isochrones are labelled 
with ages in Myr. 
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Fig. 12. — Mass-Tgff relationships from various groups (NextGen models of Baraffe et al. (1998) 
- top left; COND03 models of Baraffe et al. (2003) — top right; D'Antona & Mazzitelli (1997) 
models — bottom left; Arizona models of Burrows et al. (1997) — bottom right). The weighted 
(open circles) and unweighted (solid circles) masses of GJ 569 Ba and Bb are also plotted together 
with 1 a error bars. The vertical dotted line represents the substellar limit. Isochrones are labelled 
with ages in Myr. 
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Fig. 13. — Radius vs. mass for the models of Baraffe et al. (2003) (left panel) and Burrows et al. 
(1997) (right panel). The weighted (open circles) and unweighted (solid circles) masses of GJ 569 Ba 
and Bb are also plotted together with 1 a error bars. The radii of these objects have been obtained 
from logg = 4.75 ±0.25, which results from the synthetic spectral fitting. 
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Table 1: New astrometry of GJ 569 Bab. 



UT Date 


MJD 


Separation 


P.A. 






(") 


(deg) 


2001 Jun 28 


52088.2913 


0.1024 ±0.0012 


352.6 ±2 


2001 Sep 01 


52153.2211 


0.1033 ±0.0010 


9.7±2 



Table 2: Improved orbital parameters of GJ 569 Bab. 



Parameter This paper Previous value 3. 



Total mass, M (M ) 


0.125 


± 


0.005 


0.123 


± 


0.008 


Period, P (days) 


876 


± 


9 


892 


± 


25 


Eccentricity, e 


0.32 


± 


0.01 


0.32 


± 


0.02 


Semi-major axis, a (AU) 


0.90 


± 


0.01 


0.90 


± 


0.02 


Inclination, i (deg) 


34 


± 


2 


34 


± 


3 


Arg. periapsis, u (deg) 


257 


± 


2 


76 b 


± 


4 


Long, ascending node, O (deg) 


321.5 


± 


2.0 


141 b 


± 


4 


Epoch (MJD), T (days) 


51822 


± 


3 


51820 


± 


4 



"From Lane et al. (2001a). 

6 There was an ambiguity of 180° in the determination of these two parameters in Lane et al. (2001a). 



Table 3: Log of spectroscopic observations. 



Object 


UT Date 


Exposure 

(s) 


Airmass 


Weather 


GJ 569 Bab 


2001 May 10 


3x200 


1.00-1.01 


Clouds 


G J 569 Bab 


2001 Jun 28 


2x1200 


1.02-1.04 


Clouds 


G J 569 A 


2001 Jun 28 


2x400 


1.12-1.14 


Clouds 


G J 569 Bab 


2001 Aug 22 


3x600 


1.38-1.55 


Clouds 


G J 569 Bab 


2001 Sep 01 


2x600 


1.36-1.48 


Clear 



Table 4: Radial velocity measurements. 



MJD t>[GJ569Bb - GJ569Ba] u fc [GJ569Ba] v h [GJ569Bb] 

(kms -1 ) (kms -1 ) (kms -1 ) 

52039.452 -5.4±2.3 -9.1 ±1.1 -14.5 ±2.0 

52088.314 -4.0 ±2.3 -9.6 ±1.1 -13.5 ±2.0 
52143.272 -3.7 ±1.5 

52153.240 -2.5 ±1.5 -10.5 ±1.1 -13.0 ±1.1 
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Table 5: Best-fit masses of GJ 569 Ba and Bb. 





Vo 


M (GJ 569 Ba) 


M (GJ 569 Bb) 


Q 




(km s _1 ) 


(M Q ) 


(M Q ) 




Weighted 


-11.46 ±0.45 


0.071 ±0.011 


0.054 ±0.011 


0.76 ±0.24 


Unweighted 


-11.52 ±0.45 


0.068 ±0.011 


0.057±0.011 


0.84^ 



